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ABSTRACT 

Aims. We present an X-ray study of the supernova remnant SNR J0533-7202 in the Large Magellanic Cloud (LMC) and determine 
its physical characteristics based on its X-ray emission. 

Methods. We observed SNR J0533-7202 with XMM-Newton (background flare-filtered exposure times of 18 ks EPIC-pn and 31 
ks EPIC-MOSl, EPIC-MOS2). We produced X-ray images of the supernova remnant, performed an X-ray spectral analysis, and 
compared the results to multi-wavelength studies. 

Results. The distribution of X-ray emission is highly non-uniform, with the south-west region much brighter than the north-east. The 
detected X-ray emission is correlated with the radio emission from the remnant. We determine that this morphology is most likely due 
to the supernova remnant expanding into a non-uniform ambient medium and not an absorption effect. We estimate the remnant size 
to be 53.9 (+3.4) x 43.6 (+3.4) pc, with the major axis rotated ~ 64° east of north. We find no spectral signatures of ejecta emission 
and infer that the X-ray plasma is dominated by swept up interstellar medium. Using the spectral fit results and the Sedov self-similar 
solution, we estimate the age of SNR J0533-7202 to be ~ 17 - 27 kyr, with an initial explosion energy of (0.09 - 0.83) x 10^^ erg. 

We detected an X-ray source located near the centre of the remnant, namely XMMU J053348.2-720233. The source type could not 
be conclusively determined due to the lack of a multi-wavelength counterpart and low X-ray counts. We found that it is likely either 
a background active galactic nucleus or a low-mass X-ray binary in the LMC. 

Conclusions. We detected bright thermal X-ray emission from SNR J0533-7202 and determined that the remnant is in the Sedov 
phase of its evolution. The lack of ejecta emission prohibits us from typing the remnant with the X-ray data. Therefore, the likely 
Type la classification based on the local stellar population and star formation history reported in the literature cannot be improved 
upon. 

1. Introduction 

The study of supernova remnants (SNRs) is important in many 
fields of astrophysics. SNRs distribute the heavy elements 
formed inside their stellar progenitors and via explosive nucle¬ 
osynthesis in their supernovae (SNe) into the interstellar medium 
(ISM) of their host galaxies. The kinetic energy they impart 
drives the physical evolution of the ISM and cosmic rays can be 
accelerated in their fast moving shocks. Studies of these objects 
in the Milky Way are heavily affected by many factors includ¬ 
ing uncertain distance estimates and high foreground absorption. 

This is especially problematic for the study of evolved SNRs, 
whose low temperature plasmas are particularly susceptible to 
foreground absorption. To avoid these biases in studies of SNRs, 
we must look to the nearby galaxies of the Local Group. 

The Large Magellanic Cloud (LMC) is the quintessential 
target for the study of a population of SNRs. At a distance of 
50 kpc ( [di Benedetto ||2008] ) it is sufficiently close that its stel¬ 
lar population and diffuse structure is resolved in most wave- 

^ Based on observations obtained with XMM-Newton, an ESA sci¬ 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA 

Now at Laboratoire AIM, CEA-IREU/CNRS/Universite Paris 
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length regimes. The LMC is almost face-on (inclination angle 
of 30-40°, [van der Marel & Cioni||2001t [Nikolaev et al.||2004| ) 
making the entire population of SNRs available for study. The 
modest extinction in the line of sight (average Galactic fore¬ 
ground A^h ~ 7 X 10^^ cm“^) means optical and X-ray obser¬ 
vations of SNRs are only slightly affected by foreground absorp¬ 
tion, whereasitsloca^n in one of the coldest parts of the radio 
sky ( [Haynes et al.|199l) ) allows for improved radio observations 
without interference from Galactic emission. These favourable 
conditions have allowed a substantial populati on of SNRs to be 
identified in this galaxy. [Badenes et al.[ ( [2010[ ) compiled a cata¬ 
logue of 54 SNRs in the LMC. However, this number has been 
steadily increasing in recent years because of improved observa¬ 
tion sensitivity and/or multi-wavelength confirmation of candi¬ 
dates. Recent additions inc lude[Kavanagh et al.[([2013[) , Bozzetto 
[et al.[ ( [MT^ , [Maggi et al.[ ( |2014[ ), and [Kavanagh et al.[ ( 2015| ). 


[Bozzetto et al.[ ( [20T^ hereafter BF13) used radio observa¬ 
tions by the Australia Telescope Compact Array (ATCA) to con¬ 
firm the SNR nature of an unclassified object in the Magellanic 
Cloud Surveys of [Filipovic et al.[ ( [l995|[1998a|b| ), assigning the 
identifier SNR J0533-7202. The SNR exhibits a horseshoe mor¬ 
phology with an estimated size of 37(+l) x 28(+l) pc. This dif¬ 
fers from the dimensions determined in this work and the reasons 


for this are discussed in Section 3.1 BF13 report no significant 




































A&A proofs: manuscript no. astro-ph 



RA (J2000) RA (J2000) 


Fig. 1. Left: XMM-Newton EPIC image of SNR J0533-7202 in false colour with RGB corresponding to 0.3-0.7 keV, 0.7-1.1 keV, and 1.1- 
4.2 keV. The image is overlaid with 6 cm (5500 MHz) contours with levels corresponding to 3, 6, 9,12, and 15cr (from BF13, Fig. 1). The location 
of the point source XMMU J053348.2-720233 is shown. Right: Same as Left but with a zoomed out view showing the SNR and background 
spectral extraction regions. 


optical line emission, with only very tentative Ha emission de¬ 
tected in the Magellanic Cloud Emission Line Survey (MCELS, 
[Smith et ^|2006| ) images. Similarly, these authors find no in¬ 
frared (IR) emission associated wit h SNR J0533-7202 i n the 
Spitzer mosaic images of the EMC ([Meixner et al.||2006|). The 
ROSAT source lRXSJ053353.6-720fe4 lVoges et al.||19^ is 
located at the rim of the radio SNR. However, the source is very 
faint with only enough count statistics for a position determi¬ 
nation. BE 13 used data from the Magellanic Cloud Photometric 
Survey catalogue ( [Zaritsky et al.||2004|) and the s tar formation 
history map of the EMC ( [Harris & Zaritsky |20Q9| ) to determine 
that SNR J0533-7202 is associated with an old stellar popula¬ 
tion, free of recent star formation, and is therefore most likely 
the result of a Type la explosion. 

In this paper we report on a follow-up X-ray observation of 
SNR J0533-7202 with XMM-Newton. This observation and data 
analysis are described in Section 2. The results from the observa¬ 
tion are explained in Section 3. Einally, we summarise our work 
in Section 4. 


2. Observations and data analysis 


2.1. X-ray observation 


XMM-Newton ([Jansen et al.[[2001|) observed SNR J0533-7202 
on May 9 2013 (Obs. ID 0720440101, PI M. Sasaki). The pri¬ 
mary instrument for the observation was the European Photon 
imaging Camera (EPIC), which consists of a pn CCD ( [Striider 
et al.[[200T] ) and two MOS CCD ( [Turner et al.[[200l| ) imaging 
spectrometers. The observational data were reduced using the 
standard reduction tasks of SASQ version 13.5.0, filtering for pe¬ 
riods of high particle background. This resulted in fiare-filtered 
exposure times of ~ 18 ks for the EPIC-pn and ~ 31 ks for both 
the EPIC-MOSl and EPIC-MOS2 detectors. 


^ Science Analysis Software, see http://xmm.esac.esa.int/sas/ 
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2.2. X-ray imaging 

We produced images and exposure maps in various energy bands 
from the fiare-filtered event lists for each EPIC instrument. We 
filtered for single and double-pixel events (PATTERN < 4) from 
the EPIC-pn detector, with only single pixel events considered 
below 0.5 keV to avoid the higher detector noise contribution 
from the double-pixel events at these energies. All single to 
quadruple-pixel events (PATTERN < 12) were considered for the 
MOS detectors. 

We used three energy bands suited to the analysis of the spec¬ 
tra of SNRs. A soft band from 0.3-0.7 keV includes strong lines 
from O; a medium band from 0.7-1.1 keV comprises Ee L-shell 
lines as well as Ne Hea and Lya lines; and a hard band (1.1-4.2 
keV) which includes lines from Mg, Si, S, Ca, and Ar. 

We subtracted the detector background from the images us¬ 
ing filter-wheel-closed data (EWC). The contribution of the de¬ 
tector background to each EPIC detector was estimated from 
the count rates in the comer of the images, which were not ex¬ 
posed to the sky. We then subtracted appropriately-scaled EWC 
data from the raw images. We merged the EPIC-pn and EPIC- 
MOS images into combined EPIC images and performed adap¬ 
tive smoothing of each using an adaptive template determined 
from the combined energy band (0.3-4.2 keV) EPIC image. The 
sizes of Gaussian kernels were computed at each position in or¬ 
der to reach a signal-to-noise ratio of five, setting the minimum 
full width at half maximum of the kernels to ~ 8". In the end the 
smoothed images were divided by the corresponding vignetted 
exposure maps. Einally, we produced a three-colour image of 
SNR J0533-7202 which is shown in Eig.[T}left. 

2.3. Point sources 

We searched for point sources associated with SNR J0533-7202 
using the SAS task edetect_chain. Images were extracted 
from each of the EPIC instruments in the standard energy bands 
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0.2-0.5 keV, 0.5-1 keV, 1-2 keV, 2-4.5 keV, and 4.5-12 keV, 


with the same pattern filtering criteria outlined in Section 2.2 

For an adopted minimum-detection likelihooc0of 10, we de¬ 
tected several sources in and around SNR J0533-7202. How¬ 
ever, most of these were false detections due to the bright ex¬ 
tended emission regions of the remnant. After visual inspection 
and false source exclusion, only one was found to reside within 
the extent of SNR J0533-7202. This source is located close to 
the remnant centre, as indicated in Fig.[T]-left. 

We corrected the X-ray positions of the XMM-Newton 
sources by calculating the offsets of those stars in the field- 
of-view (FOV) with optical and infrared (IR) counterparts us- 


ing the USNQ-Bl and 2MASS catalogues of Monet et al 


( [2003| ) and |Cutri et al.| ( [2003| ), respectively. The offset between 
the X-ray positions and the optical/IR positions corrected for 
proper motion is ARA = -0.03 + 0.54 arcsec and ADec = 
0.54 + 0.80 arcsec. Since both offsets are not statistically sig¬ 
nificant, we do not apply an astrometric correction to the X- 
ray source positions. Therefore, we determined J2000 coordi¬ 
nates of RA = 05h33"^48.29^ and Dec = -72^02”^33.0^ for the 
source near the centre of SNR J0533-7202 and assign the identi¬ 
fier XMMU J053348.2-720233. The analysis and discussion of 
XMMU J053348.2-720233 and its relation to SNR J0533-7202 


is given in Section 3.3 


2.4. X-ray spectral analysis 

For the spectral analysis, we made use of the EPIC-pn and EPIC- 
MOS data. Although less sensitive, the EPIC-MOS spectral res¬ 
olution is slightly better than the EPIC-pn, and therefore helps 
to constrain the parameters of the spectral models. Before pro¬ 
ceeding with the spectral extraction we generated vignetting- 
weighted event lists for each EPIC instrument to correct for the 
effective area variation across our extended source using the S AS 
task evigweight. Spectra were extracted from elliptical regions 
encompassing the X-ray extent of the SNR. Backgrounds were 
extracted from a source and diffuse emission free region immedi¬ 
ately to the south-west (see Fig.[T]-right). Detected point sources 
in the source and background regions were excluded. All spec¬ 
tra were rebinned so that each bin contained a minimum of 30 
counts to allow the use of the;^^ statistic during spectral fitting. 
The EPIC-pn and EPIC-MOS source and background sp ectra 
were fitted simultaneously using XSPEC ( [Arnaud 1996[) ver- 
sion 12.8.1 with abundance tables set to those of I Wilms et aP 
( 2000|), and photoelectric absorption cross-s ections set to those 
of |Balucinska-Church & McCammon| ( |1992| ). 


2.4.1. X-ray background 

Detailed descriptions of the X-ray b ackground cons t ituents and 
spectral modelling can be found in |Bozzetto et al.| ( |2Q14| ) and 
|Maggi etaL] ( |2014| ). Here we briefly summarise the treatment of 
the X-ray background in the case of SNR J0533-7202. 

The X-ray background consists of the astrophysical X-ray 
background (AXB) and particle induced backgr ound. The AXB 
typically comprises four or fewer components ( [Snowden et al. 


|2008[ [Kuntz & Snowden||201Q| ), namely the unabsorbed ther¬ 


mal emission from the Local Hot Bubble, absorbed cool and hot 
thermal emission from the Galactic halo, and an absorbed power 
law representing unresolved background active galactic nuclei 


^ The detection likelihood L is defined by L = -In p, where p is the 
probability that a Poissonian fluctuation of the background is detected 
as a spurious source. 


(AGN). The spectral properties of the background AGN com¬ 
ponent were fixed to the well known values of E ~ 1.46 and 
a normalisation equivalent to 10.5 photons keV cm"^ s“^ sr“^ 
( jChen et al.|1997| ). The foreground absorbing material comprises 
both Galactic and LMC components. The foreground Galac¬ 
tic absorption component was fixed at 7.4 x 10^^ cm“^ based 
on the [Dickey & Lodamn] ( [1990[ ) HI maps, determined using 
the HEASARC Ah Tooj^ while the foreground LMC absorp¬ 
tion component, with abundances set to those of the LMC, was 
allowed to vary. We note that in all the forthcoming fits, the 
normalisation of the Local Hot Bubble component consistently 
tended to 0, with only upper limits to its contribution determined. 
Therefore, since the model is not formally required, we excluded 
this component from the spectral fits. 

The particle-induced background of the EPIC consists of the 
quiescent particle background (QPB), instrumental fluorescence 
lines, electronic read-out noise, and residual soft proton (SP) 
contamination. To determine the contribution of these compo¬ 
nents we made use of vignetting corrected EWC data. We ex¬ 
tracted EWC spectra from the same detector regions as the ob¬ 
servational source and background spectra. The EPIC-p n EWC 
spectra were fitted with the empirical model developed by [Stum 
( [2Q12[ ) and the EPIC-MOS EWC spectra with an adapted form 
of this model. Since these spectral components are not subject 
to the instrumental response, we used a diagonal response in 
XSPEC. The resulting best-fit model was included and frozen 
in the fits to the observational spectra, with only the widths and 
normalisations of the fluorescence lines allowed to vary. We also 
included a multiplicative constant to normalise the continuum to 
the observational spectra using the high energy tail {E > 4 keV) 
where the QPB component dominates. The residual SP contam¬ 
ination was fitted by a power law not convolved with the instru¬ 
mental response ( [Kuntz & Snowden|2008[ ). 


2.4.2. Source emission 

To account for the X-ray emission from SNR J0533-7202, we 
included a thermal plasma model in the spectral fits absorbed 
by foreground Galactic and LMC material. The abundance of 
the LMC a bsorption component wa s fixed to the LMC val¬ 
ues (0.5 Zq, [Russell & Dopita[[l992[ ). We made use of a non¬ 
equilibrium ionisation (NEI) model in XSPEC, appropriate for 
SNRs, namely the plane-parallel vpshock model (Borkowski 
[et al.[[200~T] ). The vpshock model features a linear distribution 
of ionisation ages behind the shock which is more realistic than 
single ionisation age mode ls such as vnei. Given its relatively 
large size (see Section [3T] ), we can assume that the remnant has 
most likely entered the Sedov phase of its evolution. Therefore, 
we also fit the spectrum with th e Sedov dynamical mod el, imple¬ 
mented as vsedov in XSPEC ( [Borkowsld et al.|2001[ ). In all of 
our spectral fits we allowed the abundances of O and Ee to vary, 
to determine if the abundances of these elements are consistent 
with the LMC ISM, or are enhanced due to an ejecta contribu¬ 
tion. The results of the fits are shown in Table [T] with the best-fit 
vsedov model shown in Eig.[^ 


3. Results 

3.1. SNR J0533-7202 morphology 

The X-ray emission from SNR J0533-7202 is highly correlated 
with the 3 cm and 6 cm radio emission detected in BE 13 with the 

^ http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/ 
w3nh.pl 
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Table 1. Spectral fit results SNR J0533-7202. See text for description 
of the models. 


Component 

Parameter 

Value 

vpshock 

phabs 

Mi, Gal (10^^ cm ^) 

0.074“ 

vphabs 

Mi,LMC (10^^ cm“^) 

0.00 (< 0.03)* 

vpshock 

kT 

0.36 (0.28-0.47)* 


0 (Z/Zo) 

0.39 (0.33-0.45) 


Fe (Z/Zo) 

0.60 (0.45-0.80) 


Tu (10^^ s cm“^) 

0.86 (0.37-2.00) 


EM (10^^ cm-3) 

4.2 (2.6-7.5) 


Fx^ (10“^^ erg s"^ cm“^) 

1.9 


Lx^(10^^ erg s“^) 

1.0 

Fit statistic 

Ay 

1.13 (312d.o.f.) 

vsedov 

phabs 

Ml,Gal (10^^ cm ^) 

0.074^ 

vphabs 

Mi,lmc (10^^ cm“^) 

0.00 (< 0.02)^ 

vsedov 

kT 

0.22 (0.17-0.36)^ 


0 (Z/Zo) 

0.48 (0.42-0.54) 


Fe (Z/Zo) 

0.67 (0.46-0.88) 


T (10^^ s cm“^) 

3.25 (0.61-25.43) 


cm-3) 

6.8 (2.7-16.2) 


Fx^ (10“^^ erg s"^ cm“^) 

1.9 


Lx^(10^^ erg s“^) 

1.0 

Fit statistic 

Ay 

1.15 (313 d.o.f.) 


Notes. The numbers in parentheses are the 90% confidence intervals 
Fix ed to the Galactic column density from the Dickey & Lockman 
( 1990 1 HI maps. Absorption and thermal component abundances 
fixed to those of the LMC. 0.3-10 keV absorbed X-ray fiux. 0.3- 
10 keV de-absorbed X-ray luminosity, adopting a distance of 50 kpc to 
the LMC. 


same horseshoe morphology and enhanced emission regions in 
the south-west (SW, see Fig.[T]-left). The X-ray emission is soft 
with no detection of extended emission in the 1.1-4.2 keV band. 
Emission in the 0.3-0.7 keV traces the outer edge of the rem¬ 
nant, except for the north-east (NE) limb, pointing to an approxi¬ 
mately elliptical shape. To estimate the size of the X-ray remnant 
we determined the average background surface brightness and 
corresponding standard deviation (cr) in the 0.3-0.7 keV band. 
We then defined the edge of the SNR as regions where the ex¬ 
tended emission surface brightness rises to 3cr above the average 
background. We then fitted an ellipse to this contour, excluding 
the absent limb emission in the NE and an enhancement due to 
a nearby point source in the SW (see Eig.|^. The error on the fit 
was determined by quantifying the standard deviation of points 
on the contour from the best-fit ellipse. We determined a best- 
fit ellipse centred on the J2000 coordinates RA = 05^33”^50.19^ 
and Dec = -72^02”^56.0^ of size 3:71 (+0:23) x 3:00 (±0:23), 
corresponding to 53.9 (+3.4) x 43.6 (+3.4) pc at the LMC dis- 
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tance, with the major axis rotated ~ 64° east of north. The best-fit 
dimensions and error are shown in Eig.[^ 

The SNR size in the NW-SE direction (minor axis) is approx¬ 
imately the same as determined in the radio analysis of BF13 
(their major axis). However, our NE-SW (major axis) size deter¬ 
mined from the X-ray data is much larger than that determined 
in BF13 (their minor axis). This is most likely due to the size 
estimation methods employed. BF13 determined the dimensions 
using a one-dimensional intensity profile along the approximate 
NE-SW axis, which does not take the elliptical shape of the outer 
edge of the SNR into account, and only considers the bright SW 
end. 

The asymmetric brightness profile in X-rays could be be¬ 
cause of either an ambient density gradient along the NE-SW 
axis, or because of a higher foreground absorbing column to¬ 
wards the NE of the remnant. To investigate this, we divided 
the remnant into two halves namely the faint NE and brighter 
SW regions, and extracted EPIC spectra. We fitted these spectra 
simultaneously in XSPEC with the thermal plasma models dis¬ 
cussed in Section [2^ For the first set of fits, we linked the fore¬ 
ground absorbing column (A^h) towards both regions and only 
allowed the normalisation of the thermal components to vary to 
determine if the SW region is intrinsically brighter. For the sec¬ 
ond set of fits we linked the thermal component parameters of 
each region, allowing only A^h to vary to determine if the bright¬ 
ness profile is due to foreground absorption. We found that, in 
both cases, the foreground absorption values were negligible and 
only upper limits to the absorption could be determined. These 
results suggest that the asymmetric brightness profile is intrinsic 
and due to an ambient density gradient rather higher foreground 
absorption towards the NE of the remnant. 

We also compared the X-ray morphology to the HI maps of 
|Kim et aL] ( [2003| ) to search for either a density gradient or vari¬ 
ations in absorption across the remnant. While the resolution of 
these HI maps (!' or 15 pc at the LMC distance) is poorer than 
XMM-Newton, they do allow for a coarse comparison with the 
SNR morphology. However, we did not find any conclusive evi¬ 
dence to support either explanation. 


3.2. SNR J0533-7202 X-ray emission 


The X-ray spectrum of SNR J0533-7202 was found to be well 
described by either the vpshock or vsedov model with little 
difference in the fit-statisitic (xl = 1.13 and;^J = 1.15, respec¬ 
tively), as shown in Table In each case, the best-fit absorption 
due to foreground material in the LMC was negligible, with only 
upper limits determined. The best-fit plasma temperatures for 
each model agree within the 90% confidence intervals, and are 
consistent with other evolved SNRs in the LMC (see |Bozzetto| 
et al.|2014{|Maggi et al.|2014"l for example). The ionisation ages 
inferred from the fits are also compatible, though the value is not 
well constrained in the vsedov model. The best-fit abundances 
of O and Fe in each model are also in agreement, and are con¬ 
sistent with the LMC abundances ( [Russell & Dopita|1992| ). This 
suggests that the thermal plasma is dominated by swept-up ISM, 
with no discernible contribution from shocked ejecta. This justi¬ 
fies our earlier assumption that the SNR is in the Sedov phase of 
its evolution, and prohibits us from suggesting the type of pro¬ 
genitor supernova responsible for SNR J0533-7202. 

Using the results of the vsedov fit we can estimate physical 
parameters of SNR J0533-7202 using the Sedov dynamical 
model (see [Sasaki et al.||2004} [Bozzetto et al.[[MTj} [Maggi 


et al. 2014[ for examples). The X-ray shell of the remnant is 


delineated by an ellipse (see Section with semi-major and 
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Fig. 2. XMM-Newton EPIC-pn spectrum of SNR J0533-7202, shown 
alone for clarity. The best-fit model is given by the solid black line. 
The red dashed line indicates the vsedov component corresponding to 
the SNR emission, the magenta dash-dot lines mark the AXB compo¬ 
nents, the blue dash-dot-dot-dot line shows the combined contributions 
of the QPB, residual SPs, instrumental fluorescence lines, and electronic 
noise. The source fit parameters are given in Table 
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Fig. 3. Combined 0.3-0.7 keV EPIC image of SNR J0533-7202. The 
red points delineate the contour level corresponding to 3cr above the 
average background surface brightness. The green solid line shows the 
best-fit ellipse to the contour, with the dashed lines indicating the Icr 
error on the fit. The white plus-sign marks the best-fit centre of the SNR 
and the white cross gives the location of XMMU J053348.2-720233. 
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semi-minor axes of 27.0 pc and 21.8 pc, respectively. Assuming 
these are the first and second semi-principal axes of an ellipsoid 
describing SNR J0533-7202, and that the third semi-principle 
axis is in the range 21.8-27.0 pc, we determined the volume (V) 
limits for the remnant and their corresponding effective radii 
(Reff) to be (1.6 - 2.0) x 10^^ cm^ and 23.4-25.1 pc, respectively. 


Table 2. Physical properties of SNR J0533-7202 derived from the Se¬ 
dov model 


no 

(10“^ cm“^) 

V, 

km s“^ 

t 

(kyr) 

M 

(Mo) 

Eo 

(10^' erg) 

2.6-7.0 

377 - 549 

17-27 

48 - 162 

0.09-0.83 


The best-fit shell X-ray temperature (see Table [T]) corre¬ 
sponds to a shock velocity 


V, 


i6^r,y 


( 1 ) 


where kTs is the postshock temperature and fi is the mean 
mass per free particle. For a fully ionised plasma with LMC 
abundances, ii = 0.6Imp, resulting in a shock velocity 

= 429 (377 - 549) km s“^ The age of the remnant can 
now be determined from the similarity solution: 



( 2 ) 


where R = Rgff and t is the age of the remnant. This gives an age 
range of ~ 17-27 kyr. The pre-shock H density (^h,o) in front 
of the blast wave can be determined from the emission measure 
{EM, see Table [T]). Evaluating the emission integral over the Se¬ 
dov solution using the approximation for the radial density dis¬ 
tribution of |Kahn| ( [l975| ) gives 


/ 


nefiY^dV = EM = 2.071 ^ 


( 3 ) 


where and riu are electron and hydrogen densities, respec¬ 
tively, and V is the volume (e.g., |Hamilton et al.|p^8^ . Tak¬ 
ing ne/nu = 1.21 and the determined range of volumes, this 
equation yields ^h,o = (2.4 - 6.4) x 10“^ cm“^. Since the pre¬ 
shock density of nuclei is given as no ~ 1.1^h,o, it follows that 
no = (2.6 - 7.0) X 10“^ cm“^. The initial explosion energy (£” 0 ) 
can be determined from the equation: 


R = 


imEof 

Enno 


. 2 \ 1/5 


(4) 


where is the mean mass per nucleus, given as = lAnip. 
This results in an initial explosion energy in the range (0.09 - 
0.83) X 10^^ erg. This is slightly lower than the canonical 
10^^ erg but consistent with several Type la SNRs in the LMC 
(see [Hendrick et al.|[200^ [Bozzetto et al.||2014| for example). 
The swept-up mass contained in the shell is given simply by 
M = V/inno, which is evaluated to (48 - 162) M©. This large 
amount of swept-up material supports our original assumption 
that SNR J0533-7202 is into the Sedov phase of its evolution. 
All the derived parameters of SNR J0533-7202 are summarised 
in Tabled 

Using the determined values of t and no we estimated the ex¬ 
pected value of the ionisation parameter (r) in the Sedov spec¬ 
tral fit to the remnant. Given that r is the product of the electron 
density immediately behind the shock front and the remnant age 
( Borkowski et al.|[2001| ), it can be written r = net ^ 4.8^h,o^- 
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We determined the expected r to be in the range (0.60 - 2.53) x 
10^^ cm“^ s, which is in agreement with the r determined in the 
spectral fits (see Table [T]), though the upper bound of the fitted 
ionisation parameter is poorly constrained. 

3.3. XMMU J053348.2-720233 

Our source detection identified the solitary X-ray 
source XMMU J053348.2-72 023 3 projected inside of 
SNR J0533-7202 (see Section ^ and Fig. [T}left). Its lo¬ 
cation near the centre of the remnant (see Fig.|^ could indicate 
that the source is a remnant compact object, pointing to a 
core-collapse origin for SNR J0533-7202. This would not 
be consistent with the Type la scenario proposed by BF13, 
and therefore the classification of this source has important 
implications for the typing of the SNR. 

We searched for a multi-wavelength counterpart in various 
catalogues in the literature but none could be identified. As such, 
the X-ray data alone were used to determine a likely classifi¬ 
cation for the source. XMMU J053348.2-720233 is embedded 
in the bright extended emission of SNR J0533-7202. We deter¬ 
mined the source and background (an annulus around the source) 
extraction regions using the SAS task region, considering only 
the 1-4.5 keV events to improve the contrast between the source 
and contaminating background emission from the SNR. Source 
spectra, background spectra, and ancillary files were then ex¬ 
tracted from each of the EPIC instruments and merged using the 
SAS task epicspeccombine. From the combined EPIC spec¬ 
trum we determined a net total of ~ 70 counts, much too low to 
perform a meaningful spectral analysis. Instead, we performed 
a hardness ratio (HR) analysis which compares the number of 
counts in certain energy bands to determine the approximate 
shape of the spectrum, allowing spectral parameters to be in¬ 
ferred. The HRs are defined as 


HRr = 


Rj+i - Rj 
^i+l + Ri 


(5) 


where R corresponds to the net source counts in a given energy 
band. For our analysis we defined the energy bands as /i = 0.3 - 
1 keV, /2 = 1 “ 2 keV, and 13 =2-4.5 keV, which provides two 
HR values. 

We determined the net source counts in the specified energy 
bands for XMMU J053348.2-720233 and calculated the hard¬ 
ness ratios with the Bayesian method described in [Park et ah] 
( |2006 j), implemented in the associated BEHR code (ver. 12-12- 
2 OI 3 Q This Bayesian method is much more accurate than the 
classical methods of HR estimation since it assumes that the 
detection of counts is a random process described by Poisson 
statistics. In addition, this method can handle the non-detection 
of counts in one or more energy bands which would cause clas¬ 
sical methods to fail. 

To determine the likely spectral parameters implied by 
the HR values of XMMU J053348.2-720233, we assumed 
that the source could be approximately described by a power 
law model absorbed by material in the Galaxy and LMC 
(phabs*vphabs*pow in XSPEC). We simulated a set of spectra 


with fixed Galactic absorption (A^H,Gai = 7.4 x 10^^ cm Dickey 


|& Loc]anan||1990| ), and varying va lues of LMC foreground ab- 


sorption (A^h,lmc; with Z = 0.5 Z©, Russell & Dopita|1992 ) and 


photon index (F). We produced a grid in HR space and plotted 
the HR values of XMMU J053348.2-720233 over this, shown 


See http: //hea-www. harvard. edu/astrostat/behr/ 
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Fig. 4. HR plot derived from a grid of simulated spectra, assuming an 
absorbed power law model, with the HR of XMMU J053348.2-720233 
shown in black. See text for a description of the HR simulation and 
chosen energy bands. The blue lines indicate models on the grid of equal 
r, whereas the red lines mark the models on the grid of equal equivalent 
hydrogen column Ah, tot = AH,Gai + Ah,lmc- The dashed red line gives 
the total Ah measured value through the Galaxy and the LMC ( Kalberla| 
|et al. 2 OO 5 I of 2 X10^^ cm"^. The dotted red line gives the co mbined Ah 
measured throug h th e Galaxy and the LMC determined by [Dickey &| 
[Lockman ( |1990| ) and Kim et al. ( |2003| , respectively, of 2.5 x 10^^ cm"^. 


in Fig.|^ The source HRs are consistent with that of an absorbed 
power law with F 1.8-3 and equivalent hydrogen column 
^^H,tot = ^H,Gai + Nu,lmc ~ (0-6 - 1.3) X 10^^ cm“^. 

In the absence of a foreground multi-wavelength counter¬ 
part, XMMU J053348.2-720233 is likely either located in the 
LMC or is a background AGN. A background AGN would 
be subject to foreground absorption by all line-of-sight mate¬ 
rial in the Galaxy and the LMC, which we estimated in two 
ways. We first obtained the value of Ah through the G alaxy and 
the LMC of ~ 2 X 10^^ cm“^ from the HI survey of 
jet al.j ( |2005| ), indicated by the red dashed line in Fig 
nativel y, summing the Galactic Ah from [Dickey & 


Kalberla 


4 Alter- 


.ockman 


1990) of 7.4 X 10^^ cm 


Kim et al. 


^ to the LMC Ah value from 
^ gives a combined absorbing column of 


(20031) 1.8 X 10^1 cm 

cm“^, indicated by the red dotted line in Fig. 
value inferred from the HR analysis sug 


- 2.5 X 10^1 
Therefore, the Ah, tot 
gests that either XMMU J053348.2-720233 is located beyond 
the LMC and is a background AGN, or is located in the LMC and 
subject to localised absorption, which could indicate an X-ray 
binary. In the absence of a high mass stellar population (BE 13), 
we would expect a low-mass X-ray binary (LMXB). As evident 
from Fig. these assertions are independent of the chosen HI 
map. From the likely spectral parameters and the count rate we 
estimate an indicative unabsorbed 0.3-10 keV fiux of the source 
to be Ax,0.3-10 keV ~ 2 X 10“^"^ erg cm^ s“\ corresponding to 
a luminosity of Lx,o. 3 -io keV ~ 6 x 10^^ erg s“^ if the source 
is located in the LMC, and could be consistent with a low lu¬ 
minosity LMXB. These values are subject to some uncertainty 
however given the low count statistics and poorly constrained 
spectral model. We attempted to perform a timing analysis on 
the source to search for X-ray pulsations, the period of which 
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could point to either an AGN or LMXB origin. However, the 
low number of net counts prohibited any conclusive results. In 
either case, it seems unlikely that the source contains a compact 
object resulting from a recent core-collapse event. 

4, Summary 

We presented the analysis of a follow-up XMM-Newton observa¬ 
tion of SNR 10533-7202 in the LMC. The main findings of our 
analysis can be summarised as follows: 

- We detected bright X-ray emission from SNR J0533-7202. 
We estimated the size of the remnant to be 53.9 (±3.4) x 
43.6 (±3.4) pc, with the major axis rotated ~ 64° east of 
north. The X-ray emission is highly correlated with the radio 
emission detected by BF13, displaying a highly asymmetric 
brightness distribution. We found that this is most likely due 
to an ambient interstellar medium density gradient. 

- The X-ray spectrum of the remnant is thermal in origin 
and well described by both the non-equilibrium ionisation 
vpshock model and the Sedov dynamical model. The abun¬ 
dance patterns determined from the spectral fits are con¬ 
sistent with swept-up interstellar medium and show no en¬ 
hancements due to shocked ejecta, consistent with a remnant 
in the Sedov phase. The plasma temperature is relatively soft 
(~ 0.2 keV for the Sedov model), similar to other evolved 
remnants in the LMC. 

- Using our spectral fit results and the Sedov self-similar so¬ 
lution, we estimated the age of SNR J0533-7202 to be 
~ 17-27 kyr, with an initial explosion energy of (0.09 - 
0.83) X 10^^ erg, and a swept-up interstellar medium mass of 
48-162 Mo, further confirming that the SNR is in the Sedov 
phase. 

- We determined that XMMU J053348.2-720233, the source 
located near the centre of SNR J0533-7202, can either be 
a background AGN or a source in the LMC, most likely a 
LMXB. Because no shocked ejecta emission was detected in 
the X-ray spectrum, we cannot offer a typing of the progeni¬ 
tor supernova, and the assertion of BF13 that the remnant is 
the result of a Type la event cannot be improved upon. 
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